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Abstract 

The work is devoted to a consideration of the different determinations of the proton charge radius. It is 
demonstrated that the results from the elastic electron-proton scattering have to be of a higher uncertainty. 
A review of the hydrogen Lamb shift measurements and the radius determination from them is also presented. 

1 Introduction 

The proton is the hghtest and simplest stable hadronic system. Investigation of its structure is quite important. 
The charge radius determined by the charge distribution inside the proton is one of the universal fundamental 
physical constants, because it is important for a number of very different physical problems. Some recent 
precise results are collected in Table ^ and also presented in Fig. 1. The proton charge radius is defined here as 

= \/ (^^)- Oiis can see that the values are obtained from different branches of physics. We give first a short 
description of all well known and most 'popular' results summarized thereQ 

The references and methods used to determine the proton radius are: 



*E-maiI: sgk@onti.vmim.spb.su 
t Temporary address 
-^E-mail: sek@mpq.mpg.de 

^Of course, the Lamb shift result mentioned there is too fresh to be 'popular' but it is going to be. 



Value 


Reference 


Method 


0.809(11) fm 
0.862(12) fm 
0.64(8) fm 
0.88(3) fm 
0.847(9) fm 
0.890(14) fm 


Stanford, 1963 
Mainz, 1980 
Draper, Woloshin, Liu, 1990 
Leinweber, Cohen, 1993 
Mainz, 1996 
Garching, 1997 


scattering experiment & empirical fitting 
scattering experiment & empirical fitting 

lattice QCD in chiral limit 
lattice QCD & chiral perturbation theory 
dispersion relation fitting 
hydrogen Lamb shift measurements 



Table 1: Proton charge radius 
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Figure 1: The proton charge radius as it is usuaUy known 



• In paper [y resuhs on electron-proton elastic scattering were obtained and fitted using a simple empirical 
formula. The experiment has been done by means of the accelerators and particle physics. Theoreti- 
cally they are based on the particle physics phenomenology and the quantum electrodynamics theory for 
corrections. 

• In work j2j new results for e—p were presented. The radius was found by an extrapolation to zero transfer 
momentum. 

• In paper ||^ (see also Ref. Q) a lattice calculation was performed within the chiral limit, where the pion 
is massless (m^ = 0). The result for the radius was found by fitting of the form factor computed there. 
The lattice calculation are result of the quantum field theory without using of the perturbative expansion 
in the Euclidian space. As far as the space is not continuous there are a lot of problems there due to a 
formal description of relations between the physics in the infinite continuous Minkowsky space and in the 
discrete finite Euclidian one. 

• In article |^ a chiral perturbation correction to previous calculations was presented. The chiral perturba- 
tion theory is another branch of the quantum field theory where the smallness of the masses of the u and 
d quarks is used. 

• In Ref. a many parameters fit was performed. It included transferred momentum within a wide 
range. At least two new experimental subfields are involved into this approach. First the wider range 
of data of the proton electron scattering includes a lot of absolutely different experiments. And second 
the dispersion approach based on the theory of the analytic properties of the scattering amplitude and 
other values involves a number of relation between the data from other kinds of collision. For instance, 
the neutron data obtained from scattering on nuclei are included in the evaluation. 

• The measurement of the Is — 2s transition frequency in the hydrogen atom |^ gives a result for the 
ground state Lamb shift and hence for the proton radius. The precise investigation of the hydrogen 
spectra involves laser spectroscopy to measure the transition frequencies and quantum mechanics and 
quantum electrodynamics to calculate all contributions to the energy levels beside the radius term. 

These short explanations of the results demonstrate the universality of the proton charge radius. This gives 
us the possibility to perform some cross-checking of a number of different pieces of our knowledge on the atoms 
and particles. There are also a number of other subfields where the knowledge of the radius is interesting, but 
they cannot give any precise result. An example of those can be some models of the proton or neutron. As 
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some important applications the proton polarizability and the empiric determination of the neutron magnetic 
radius can be mentioned. 

The results are given in Table 1 and Fig.l within the chronological order, but we will consider them in 
detail starting with the Lamb shift measurement. Next we will discuss the lattice calculations and then the 
scattering data and their different fitting. As one of the applications the hydrogen hyperfine structure and so 
called polarizability contribution is considered in Appendix. 

2 The hydrogen Lamb shift for the proton radius 

The nuclear dependent correction to the Lamb shift in the hydrogen atom as it is well known is of the form^ 

O Tl 

If the energy were measured with high accuracy and all other contributions to energy were known one could 
extract the proton charge radius from hydrogen spectroscopic data. However, it is only possible to measure some 
differences of the state energies: either splittings or transition frequencies. They are dimensional values and to 
determine them one has to use some other dimensional reference values. We give a review of all experimental 
results, but before that we would like to attract the attention to a problem appeared recently due to a significant 
progress in the optical measurements. 

There are two kinds of the optical measurements from which one can extract the value of the Lamb shift. 
First let us consider absolute optical measurements. They have been done for the Is — 2s, 2s — 8s and 2s — 8d 
transitions. One problem to utilize their result is unknowledge of the Rydberg constant which mainly determines 
the value of the transition frequency. The Lamb shift may be found only by using experimental results for two 
or more transition and solving a system of equations where the constant is one of the variable to find. In this 
way one gets the value of the Lamb shift from some special difference of two frequencies where the Rydberg 
constant contribution is canceled: 

^ - ^) • ^(1^ - 2^) - - ^) • ^(2^ - 8^/^) • (1) 

The other method of optical measurements (so called 'ground state Lamb shift measurements') is based 
on producing difference like in eq.(|l]) experimentally as a beat frequency. Measurements of Is — 2s/2s — 4s, 
Is — 2s/2s — 4p, Is — 2s/2s — M, Is — 3s/2s — 6s and Is — 3s/2s — Qd transitions have been done. 

One can see that in the optical measurements it is not possible to measure the Lamb shift of one separated 
level. Usually they can determine something like 

AEl{1s)-CAEl{2s) + ... , 

where the known constant C depends on the experiment and usually it lies between 4 and 5. The points stands 
for the Lamb shift of higher s levels and of levels with higher orbital momentum I > 0. All s level energy depend 
on the proton radius. 

A way of evaluating the data in which a combination of Lamb shifts of different s-states is involved has been 
considered in our works |^. It is possible to recalculate all ns shifts to the ground state value by using an 
auxiliary difference 

A(n) = AEl{Is) - n^AEL{ns) , 

which is radius-independent. The details of this calculation can be found in our paper pO[ |. 

The higher I levels are not radius-dependent on this level of accuracy and they can be also calculated with 
good enough precision ^. 

2.1 Status of the Lamb shift theory 

We now can compare the status of the ground state Lamb shift theory and the A{n) theory. All recent results 
for the Is (or 2s) Lamb shift are presented in Table 2. In the Table we have included only results obtained 
directly for the hydrogen atom. However, it has to be mentioned that some works for hydrogen-like ions with 
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Order 


Reference 


a{Za)^m 


Pachucki, 1993, ||ll| 


aiZaYra log Za 


Karshenboim, 1994, g 


a^{Zafm\og^ Za 


Karshenboim, 1993, 


a'^{Za.)^m 


Pachucki, 1994, || 




Eides, Shelyuto, Crotch, 1995, 


[Zafm^ iMlogZa 


Doncheski, Grotch, Erickson, 1991, [p^ 




Khriplovich, Milshtein, Yelkhovsky, 1992, [l|] 




Fell, Khriplovich, Milshtein, Yelkhovsky,1993, @ 




Pachucki, Grotch, 1995, |l|] 


{Zafm^/M 


Pachucki, Grotch, 1995, ^ 




Yelkhovsky, 1996, ^ 




Eides, Grotch, 1997, [|o| 




Yelkhovsky, 1997, |l| 


a[Zafm'^/M 


Bhatt, Grotch, 1987, [|2| 




Pachucki, 1995, || 




Eides, Grotch, 1995, [|| 



Table 2: New results for the Is and 2s Lamb shift 
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Order 


Reference 


Contribution 


[Zafra^lM 


Pachucki, Grotch, 1995 
Artemyev, Shabaev, Yerokhin, 1995 
Yelkhovsky, 1996 
Eides, Grotch, 1997 
Yelkhovsky, 1997 


-7.4 kHz 
-7.1(9) kHz 
2.8 kHz 
-7.4 kHz 
-16.5 kHz 


a{Zaffn'^/M 


Bhatt, Grotch, 1987 

Pachucki, 1995 
Eides, Grotch, 1995 


-20.3 kHz 
-13.9 kHz 
-20.7 kHz 



Table 3: Theoretical discrepancies for the Is Lamb shift 



the nuclear charge Z can be useful also for Z = 1. Particularly, some important numerical results were obtained 
in works [|| (am, all order of (Za), Z = 5,10...) and H {m'^/M, all order of (Za), Z = 1,2...). 

One can see from Table 2 that some corrections were investigated by a number of authors and it has to be 
mentioned that the results were not always in agreement. For instance it was a long and dramatic history of 
the calculation of the logarithmic corrections {Za)^m'^ /AI log Z a which after all were found to be equal to zero 
i pTf . Two actual contradictions of the ground state Lamb shift are presented in Table 3. For the interpretation 
of the numerical results of Ref. p^ ], which include also all higher order corrections, it is necessary to remember 
that uncertainty from the higher order corrections can be estimated as the {Zafrr? jM term multiplied by 
(Za)^ logZa, or on level of 1 kHz. 

Now we can consider the status of the specific difference A(n) = AEl{Is) — n^AEL{ns) and the Lamb 
shift of the p states. One can see from Table 4, that the theoretical status of A(n) and p state energy levels 
is physically more clear than the status of the ground state Lamb shift. Some details can be found in Refs. 

The general expression of the difference A(n) is of the form |jlffl 

A(n) = ^^(^4 X I - 1 log (i + 



(ZaV X 



3 koins) V M 
(^4( logn - ^{n + 1) + v^(2)) - log 



14 Zm 



ip{n + 1) - tp{2) - logn -I- 




where \ogko{ns) is the Bethe logarithm and V'(^) = (d/dz) logr(z). 
Here: 

• G^^{Za) is the one-loop self-energy correction of order a{Za)^Tn and higher; 

• AQff(n) is the a(ZQ;)®m-contribution of the vacuum polarization | |29[ |; 

• Bq2 is the leading logarithmic two-loop correction coefficient |3^, |2^ . 

The final results for the difference A(n) in the hydrogen and deuterium atoms as well as for their isotopic 
shift were found in Ref. p0| and they are summarized in Table 5. 

^We use a convention: h = c = 1, a = . The nuclear charge Z is useful to classify the QED contributions and in the hydrogen 
atom indeed one has Z = 1. We always give expressions for E as for the energy, but the numerical values (in kHz) are given for 
the corresponding frequency value f = E/h = E/2it. 
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Order 




A(n) 


EL{np) 




coefficient w 30 


coefficient w 1 


coefficient w -1 


a{Zafmlog^ Za 


unknown 








a^{Zafmlog^ Za 


30 kHz 








a~(Za)^' >u lop,'" Z(\ 


UlllvllOWll 


known 


known 


a^{ZaYm 


unknown 





known 


[Zafm^/M 


disagreement 





agreement 


a{Zafm^lM 


disagreement 








{Zafm{mRj,f 


radius-dependence 









Table 4: Comparison of status of -El (Is) and A(n) 



n 




A^^"(n) 


A^"°(n) 




[kHz] 


[kHz] 


[kHz] 


2 


-187232(5.5) 


-187225(5.5) 


7.3 


3 


-235079(10) 


-235073(10) 


5.9 


4 


-254428(12) 


-254423(12) 


4.7 


5 


-264162(15) 


-264158(15) 


4.0 


6 


-269747(15) 


-269743(15) 


3.5 


7 


-273246(16) 


-273243(16) 


3.2 


8 


-275583(17) 


-275580(17) 


3.0 


9 


-277221(18) 


-277218(18) 


2.9 


10 


-278413(19) 


-278410(19) 


2.7 


11 


-279308(19) 


-279305(19) 


2.6 


12 


-279996(20) 


-269993(20) 


2.5 



Table 5: Results of A(n) for n = 2 ... 12 
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By using the differences considered above one can easily recalculate the Lamb shift of the Is state to the 
Lamb shift of the 2s. In this work we prefer to discuss the Lamb splitting between the 2sx/2 smd 2pi/2 states. 
Our motivation to do that is the possibility to compare the data evaluated with A(n) (the optical measurement) 
to the one which obtained directly (the radiofrequency measurement). 

Above we have compared the status of the 2s Lamb shift and the difference A(n) and given the complete 
theoretical expression for the second jl^. The theory of the 2s Lamb shift is much more complicated. We 
mainly agree with a consideration in review |3ll . However, it has to be mentioned that compared to it, our 
result is shifted by -3.6 kHz because the a^(ZQ;)®mlog'^ Za-correction has not been included there. We 
would also like to explicitly present here the sources of the theoretical uncertainty: 

• unknown a{Zaym and higher order corrections are estimated to 1 kHz; 

• a^{Za)^m\og'^ Za and higher order terms can give up to 2 kHz; 

• the o?(Z(xy'm contributions previously estimated here to 2 kHz need more understanding. 



2.2 Status of the Lamb shift experiment 
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Figure 2: Different values of the Lamb splitting i?(2si/2) — E{2pi/2) 

All important result for the Lamb splitting (2si/2 ^ '^Pi/2) are presented in Fig. 2. It is helpful first to 
describe all labels in the text and next to discuss them. The references (and the transitions) of the experiments 
and theory presented in the picture are: 

Opt. Beat. Fr. Optical beat frequency measurements were realized with a high accuracy only recently: 

G Measurement of Is - 2s/2s - 4s, Garching, 1994: H. 
Y Measurement of Is - 2s/2s - 4p, Yale, 1993: |||. 
P Measurement of Is - 2s/2s - 4p, Paris, 1996: pi. 
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FS Fine structure measurements deal with tlie radiofrequency splitting (2^3/2 — 2si/2) 



HP The best and most recent result is of E. W. Hagley and F. M. Pipkin, J 994 |35|. 

Old Older experiments were not very precise and it is more convenient to consider them altogether. They 
include results obtained by 

o B. L. Cosens and T. V. Vorburger, 1969 ^; 
o T. W. Shyn et al, 1971 ^\ 

o K. A. SaGnya et al., 1980 @ (result corrected according to E. W. Hagley and F. M. Pipkin, 



1994 |35 ) 



LS Lamb splitting measurements have been performed directly for the classical splitting 2si/2 — 2pi/2- 

LP The result by S. R. Lundeen and F. M. Pipkin, 1981 is one of two relatively recent ones. 

SY/K The other experiment was performed by Sokolov and Yakovlev and the result corrected by us. The 
experiment properly by Yu. L. Sokolov and V. P. Yakovlev, 1982 gives a result for the ratio 
of the 2pi/2 decay rate and the splitting. The correction to decay rate was found in Ref. H, ^ 
by S. G. Karshenboim, 1994. See also Comments to On the Accuracy of Lamb Shift Measurements 
in Hydrogen, (Physica Scripta, 55 (1997) 33-40) by V. G. PaVchikov, Yu. L. Sokolov, and V. P. 
Yakovlev ||. 

Old Older experiments are summarized in one point. The results taken into account include values 
obtained by: 

o D. A. Andrews and G. Newton, 1976 pSj; 
o R. T. Robiscoe and T. W. Shyn, 1970\^. 

Theo Some theoretical values are presented too. Different results appear due to the different values of the 
proton radius discussed above. They are labeled with the proper radius values (see Table 

.862 - Mainz scattering result 
.805 - Stanford scattering result 
.847 - Mainz dispersion fitting result 

Abs.Fr.(G1997) The comparison of absolute optical frequencies is presented separately from other experimental 
points because in contrast to them it contains data from two independent experiments. The other reason 
is that this gives now the results with the highest accuracy (except the Sokolov- Yakovlev result which is 
discussed below). The Is — 2s measurement Q done in Garching (1997) cannot be used alone. The other 
absolute measured value used to find the Lamb shift is a Paris result for the 2s — 8s/8d frequency. We 
give two points due to two different Paris experiment: 

P1993 Earlier result: Paris-1993 
P1997 Recent result: Paris-1997 [|||. 

Aver It is also convenient to compare results of Abs.Fr.(G1997) with some average values. We give average 
values for three kind of experiments. The labels are easily to recognize. We would like to comment them: 



OBF The evaluation is performed with by A(n) |J( 
FS The result is obtained without using A(n), but with applying the theoretical results for the fine 
structure interval 2pi/2 — 2^3/2 J2^/ . 

LS The result is obtained without using either A{n) or the fine structure 2pi/2 — 2^3/2- The SY/K- value 
is not taken into account here. 

Gr.Av. The grand average value is found over the three average values mentioned above. 

Before a detailed discussion we would like to attract the attention to an important point. As it is 
mentioned the difference A(n) has been involved into evaluation of optical beat frequency {Opt. Beat. Fr) and 
the absolute frequency (Abs.Fr.) measurements the data. The fine structure (FS) and the Lamb splitting (LS) 
measurements data are needed no results of A(n). One can see in Fig. 2 that these two group of the results 
(A(n)-dependent and A(n)-independent) are in fair agreement. 

Now we are listing some problems for the hydrogen measurements adjustment needed to be solved for 
extracting the proton radius from the Lamb shift data. 
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The first problem to solve is to fix the experiments to be included and the uncertainties to be used for 
them. 

o Sokolov and Yakovlev experiment [^0[ and its uncertainty is an open question because some criticism 
of Hinds [|7|; 

o Older experiments altogether can give an uncertainty of about 10 kHz. Not all of them are taken 
into account in Fig. 2. We would like here to mention especially some Lamb splitting experiments 
and one for the fine structure. 

— The value of the Lamb splitting |Q is included and the uncertainty {S — 42 kHz) is taken by 
us according to the original work in contrast to a consideration by B. N. Taylor et al, 1969 
(the uncertainty there is 5 = 64 kHz). 

— The result of the Lamb splitting by S. Triebwasser et al., 1953 is excluded because the 
work includes results for two lines. Being corrected (see [Q for details) they are inconsistent; 
d = 64 kHz (according B. N. Taylor et al. @). 

— The experiment by Yu. L. Sokolov, 1973 for the Lamb splitting is not included being in 
straight disagreement with all other measurements {6 = 40 kHz). 

— The result for the fine structure by S. L. Kaufman et al., 1971 |^ is excluded also being in 
straight disagreement with all other measurements {6 — 40 kHz) . 

One can see from the picture that the two absolute measurement of the 2s — 8d/8s transition from Paris 
are in a disagreement: the earlier result of 1993 is 

/(2s - 84/2) = 770 649 561567(10) kHz 

and the recent result of 1997: Q is 

/(2s - 84/2) = 770 649 561 585(5) kHz. 

Actually these are two different experiment where different standards have been used and the results are 
expected to be partly independent. In another work the older experiment result is presented as 

/(2s - 84/2) = 770 649 561 571(12) kHz. 

It is not clear if we have to use one of them (which one ?) or to find some average value. 

It has to be also mentioned that actually there are a number of older less precise measurements from Paris 
which are quite correlated. 

Another question that appears for such high accuracy as 3 kHz for the 2s Lamb shift is the fine structure 
constant to be used. The different choice of it can lead to shift of the £'(2si/2) — E{2pi/2) splitting value 
up to 1 kHz for the results obtained by the optical beat frequency and absolute frequency measurements. 
For the £'(2si/2) — E{2pi/2) value from the fine structure measurements the shift is up to -2 kHz, where 
'-' indicates that the shift is in the opposite direction to the one for the results extracted from the 
Opt.Beat.Fr. and Abs.Fr. experiments. The results for the direct measurements of Lamb splitting cannot 
be significantly shifted. That is because all changes due to correcting of the fine structure constant value 
come via the (Za)'^ m term (the relativistic correction to the Schrodinger energy which is also responsible 
for the hydrogen fine structure). This term is about an order of magnitude larger than the Lamb shift 
properly. Due to these different signs the grand average value is less sensitive to the choice of a, but 
according to this choice the set of results may be more or less inconsistent. 

One can combine the data and extract the value of the fine structure constant. The uncertainty is of the 
level of 3 • 10^^ and it seems that it cannot be reduced in the nearest future. 

ome new problem can appear when taking into account the deuterium results. It is important to remember: 

Isotopic shift of the Is — 2s transition frequency has been measured in Garching within 0.15 kHz uncer- 
tainty n. 

Isotopic shift of A(n) with uncertainty less than 0.05 kHz has been obtained by S. G. Karshenboim, 1997 
[0 (see Table I). 
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• That gives the isotopic shift 



o of the IS Lamb shift with an uncertainty of 0.17 kHz; 
o of the 25' Lamb shift with an uncertainty of 0.02 kHz. 

• The hst of the deuterium results similar to the case of the hydrogen has to include a number of old and 
recent ones. They are shortly discussed below^. 

o The measurement of the Lamb splitting by B. L. Cosens, 1968 
hydrogen one {S = 33 kHz). 

— The experiment by S. Triebwasser et al, 1953 for the Lamb splitting also has contained correla- 
tions with the hydrogen result and leaded to the same problem as for hydrogen: the results for two 
lines corrected according to B. N. Taylor et al, 1969 |Q are inconsistent {S = 35 kHz). 

— The work by E. S. Dayhoff et al, 1953 |^ devoted to the fine structure have contained a disagreement 
between the results for two lines {6 = 53 kHz). 

o The optical beat frequency measurement of Is — 2s/2s — 4s transitions in the deuterium atom (Garch- 
ing, 1994, has been done with the same experimental setup as for the hydrogen {6 = 28 kHz). 

o The best absolute measurement has been performed recently for the 2s — 8s/8d transitions in Ref. 
p6[ (Paris, 1997) and it has been correlated with the hydrogen measurement; the uncertainty is the 
same as for the hydrogen {S — 5 kHz). 

• One can see that all deuterium results are more or less correlated with the hydrogen results obtained by 
the same teams. 

• It seems that the Sokolov and Yakovlev experiment is the only direct way where it could be possible to 
obtain uncertainty below 1 kHz for the 2s Lamb shift. Due to a perfect reproductibility, the unknown 
systematic error can be expected to be the same for hydrogen and deuterium. 



54 has been correlated with the 



3 Lattice calculations 

One of the other methods that gives a value for the proton radius (see Table |]and Fig. 1) is the approach based 
on the lattice calculation within the chiral limit in which the pion mass is equal to zero. We expect that the 
result of Ref. Q included in the Table 1 and Fig. 1 is not quite correct. Our opinion is based on the following 
critical remarks: 

• Small lattice (24 x 12 x 12 x 24) has been used in the calculation. Some progress is possible: in a recent 
evaluation of the neutron form factor in Ref. a lattice 20'^ x 32 was used. In case of the small lattice 
only a not too large number of discrete momenta are available. The proton has to occupy a relatively 
large portion of the volume of the whole space. 

• That is not an ab initio calculation. It includes a number of free auxiliary parameters found by a barions 
masses fitting. 

• The chiral limit (ttitt ~ 0) also leads to some limitation for the available values of the low momentum. 

• The finite space implies neither = nor — > 0, because only some discrete value of q are allowed. The 
radius cannot be found directly, but only by some fitting. 

• As far as the form factor Gsiq^) could be found only at not too small q^ (e.g. 0.16 Gev^ ~ Ifm^'^), it 
has been fitted with monopole or dipole fits. They have obtained results for only a few values of and 
used a normalization at = 0. 

• The result looks very sensitive to chiral perturbating: the value of the proton radius was shifted from 
0.64(8) fm [§ to 0.88(3) fm |]. It may seem quite strange that using the chiral correction a smaller 
uncertainty has been obtained. 

• But a more important question is the status of their estimation of the uncertainty. 

^The items marked by ' — ' should rather be excluded from adjustment in contrast to ones with 'o' which should be taken into 
account. 
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o Due to the simple Yukawa picture for the proton with the charged pion clouds it is not quite clear 
how the proton charge radius can be finite in the massless pion limit. The virtual neutral pions also 
lead to a divergency which is similar to the well known divergency in the electric form factor of the 
electron in the usual quantum electrodynamics. 

o The proton radius is also expected to be infinite in the pure chiral limit according to the correcting 
fit of Ref. 1). 

o But in the previous work ||^ the same authors just used only the chiral limit and obtained a finite 
result ! 

o If they are not able to recognize an actual divergency how one can believe their result and especially 
the uncertainty ? 

o We expect that was possible maybe because they used only a few of the values of not too low momenta 
and the dipole or monopole fits. The divergency were actually cut off with the finite volume of 
the space. It also has to be mentioned that the normalization G{Q) = 1 was also applied by them, 
but it has no clear physical meaning because in finite space ^ 0. 

o The results are actually sensitive to the fitting procedure, which can allow to obtain a finite value 
instead of infinity. 



4 Scattering data 



We now can start to discuss the other traditional way to investigate the proton radius based on the unpolarized 
elastic electron-proton scattering. The cross-section described by the well known Rosenbluth formula is of the 
most simple form in terms of the Sachs form factors: 



da{Eo,t 



da{Eo, 



(2) 



where 



daiEo,i 
dfl 



9 ? 

a cos 



4:E^ sin^ § 1 + ^ sin- 



2 9 
2 



is the cross-section for scattering with a spinless point-like particle (instead of the proton) and q^ — qo—^l^ < 0. 
The coefficients in eq.(^ are known 

A{q') = ' 



and 



B{q',9) = 



1 -q^/AM^ 
1 



1 -g2/4Af2 



2 tan^ 



In the limit of the small momentum q^ they are equal to 

A(0) = 1, 

5(0, 9)^1 + 2 tan^ 



The proton (electric) squared radius in terms of the form factor is defined as 

1 dGE{q^) 



n2 

P QGe{Q) dq^ 



(3) 



The Sachs form factors Ge and Gm are not only useful in the particle physics. The others for particle with 
spin 1/2 are the Dirac (Fi) and the Pauli (F2) form factors. In their terms the vertex of the emission of a 
virtual photon by a real proton is of the form 



T^ = Fi{q')-/^ + F2{q')l-C7^, 



2M 



11 



There are two relations between the different form factors: 



^'^'^ > ~ l-gV4AP 

and 

p^^q2^^_GE{q')-GMiq^) 



By measuring the scattering cross-section and applying the Rosenbluth formula one can obtain the values 
of the form factors. A number of problems appear when fitting the scattering data. 

• Dipole approximation. 

The so called dipole approximation is of the form 



GE{q') = ^^ 
1 + K 



A2 



A^-q 



This is a very good approximation for any preliminary estimation, but it has no theoretical status. It was 
used in some old fitting. 

Extrapolation . 

In the Mainz experiment 1^ they used a three term extrapolation 

GE[q ) = — ;— = ao + oi r + ^2 9 • 

1 + K 

The same was done in some older works too (see e. g. [^l[). From the charge normalization one could 
expect oo = 1. But actually the problem is not so simple and we discuss that below after consideraing 
the problem due to the magnetic form factor. 

Magnetic radius. 

In the Mainz- 1980 evaluation of the data as well as in a lot of previous works it a relation 

GE{q') = 

1 + K 

was used for fitting. That definitely works for = 

GEiO) = ^ = 1. (4) 

and definitely does not work in general because, e. g. 

G£;(4M2) = Gm(4M2). 

Hence, this is only an approximation and it is not quite clear what systematic errors come from the 
difference between GE{q^) and GAf(9^)/(l + k). In the Mainz experimental work |^ they used for the 
evaluation only data for which the magnetic effects are less than 10%. That means that the uncertainty for 
Rm cannot be smaller than 10 times 5 Re, if we have nothing else to estimate the difference Re — Rm- In 
the Stanford work they used also some higher momentum measurements. We give all accurate experimental 
points for the electric form factor in two pictures. In Fig. 3b results are included for range of momentum 
used in the Stanford fitting (up to = 3fm~^), but in Fig. 3a we keep only the same momenta as in the 
Mainz experimental fitting for which magnetic effects can be neglected (q^ < 1.4 fm~^). Experimental 
status of the problem can be illustrated with a result by H. Treissen and W. Schiitz, 1974 |^ done 
for the not too high momenta q^ 

GE{q') = 1.0m.^^. (5) 

In case of a low momentum transfer one can write 

GE{q')»l + ^ (6) 
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and 



GM(g2)« (1 + k) (^1 + 



6 

and for all low energy of the Mainz data that they used for the fitting it was 

1.5- 10-2 < (ge-i\ « < 15 • 10-2 . 



(7) 



That means that for so low momenta the result of eq.(|^) comes almost from the normalization eq.(^ 
and actually cannot give any information on the radii. 

For high there are no direct connection between the radii and form factors, but for the low momenta 
as we can see the electric radius term is of the order IQ-^. Thus there is no significantly accurate result 
for critical comparison of Re and Rm ■ 

• Normalization of cross section. 

Measuring cross section one has to take care of a proper normalization factor. Due to the experimental 
normalization the result of extrapolation 

GE{q'^) = ao + aiq^+a2q^. (8) 

to has not to give oq equal to 1 exactly. That is not possible due to the measurement nature of the 
normalization. But of course oq has to be in agreement with 1 within the experimental uncertainty. 
And an more point: the uncertainty mentioned has to be interpreted as the systematic error and thus the 
value of uq can be different for different experiments. E. g. in the Mainz-1980 work they used different 
values to fit the data of three experiments: 

- ao{Mainz) = 1.0014 [g; 

- aoiOrsay) = 1.020 @; 

- ao{Saskatoon) = 1.008 |6l|| . 

The physical meaning of the constant oq is: the function G{q^) determined straightforwardly from the 
Rosenbluth formula is not the true form factor. The true one is G{q^)/aQ. The true value of radius has 
to be defined as 

Rp - 7. — • 

From the Mainz experiment article |^] it is not quite clear if they used the definition above or 

The last equation is wrong. That can be important for incorporation of the Orsay and Saskatoon data. 
In the dispersion fitting work Q they have not used different normalizations for different experiments, 
but only 

Ge{Q) = 1. 

As a result they have underestimated the uncertainty because all systematic errors are divided there by a 
large statistical factor due to the number of the experimental points. 

We would like to discuss here shortly a well known contradiction between the Stanford ||l| and Mainz 
empirical fitting results for the radius. The Mainz result was obtained by evaluating Mainz data properly B], 
Saskatoon data and Orsay data |6^. The Stanford radius was found by treating the Stanford |]l|, Orsay 
p0| and some less precise data. The last are not presented in Figs. 3. One can easily see from the Fig. 3b, that 
the data used in the Mainz work are much anymore representative and significantly more precise. Some high 
momentum points in Fig. 3a cannot really change the situation because they are not too accurate. We will not 
consider more the Stanford data here. 

An other important point is due to the QED corrections. In order to determine the true proton form 
factor from the Rosenbluth formula one has to take into account a number of such corrections: 
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• The main radiative corrections are due to the electron form factor and the bremsstralung, and the electronic 
vacuum polarization. The next important term is the two photon exchange. The complete result in the 
one loop approximation was found in Ref. |]6^ . However, in earlier experimental works some incomplete 
results (see Refs. |6^, ^) were used. E. g. the QED corrections according to Ref. were used in 
papers but in the evaluation in work the results of Ref. |^ were used. Before including the 
data mentioned above to any compilation it is necessary to correct them according to Ref. [|64[ . It is not 
clear if that has been done or not in the both Mainz compilations. 

• Two-loop radiative corrections of the order o? include large squared logarithms log(g/me). In Ref. 
which is now used as the most complete QED result the authors claimed that part of the two-loop 
corrections may be taken into account using an exponential of the one-loop correcting factor b\ 

1 + (5i ^ exp(5i) . 

That may allow to include some leading logarithmic terms. In earlier experimental works they were not 
taken into consideration and it is not clear if they have been actually included in the Mainz experiment 
evaluation g. In the Mainz-1980 work there is no prescription how the QED corrections have been 
treated. The two-loop correction may contribute on the 1% level of the {G — l) value for a few lowest 
momenta of the Mainz experimental data. The result for the two-loop vacuum polarization is well 
known and the two- loop form factor of the electron was found in Ref. . 

• The muonic and hadronic vacuum polarization effects can be also important. The muonic and hadronic 
vacuum polarization leads to a correction lying in a range from 1% to 0.5% of the G — 1 value for 
all experimental points in Ref. 0. It has to be compared with about the 1% precision in the radius 
determination. 

• In Mainz-1980 there are no explicit equations for the Rosenbluth scattering. In earlier works often the 
massless electron approximation mg — was used. The correction is of the order ml^/ E"^ and it has to be 
lO""* — 10~^. It is not clear if older data are corrected and if in Mainz-1980 data this approximation has 
been used or not. Being £'o-dependent it cannot be incorporated in the fit. 

Thus our QED summary contains two statement 

• It is not clear if the older data from Refs. ||6^, ^ are included in the final Mainz compilations being 
corrected according to Ref. |Q or not. 

• All QED corrections which can be important on the 1% level precision of the charge proton radius deter- 
mination have been known and they have to be taken into account when fitting the data. 

Now we would like to discuss some results from the dispersion fitting approach. It allows to incorporate 
data with both high and low momentum transfer. The results for electric and the magnetic proton radii as well 
as for the magnetic neutron radius have uncertainties on the level of 1% or 0.009 fm Q. That seems strange 
because: 

• The low energy fitting and high energy one can be done absolutely separately. 

• It is clear that we have some additional information only for the proton electric radius which comes 
from the low energy scattering when the magnetic form factor effects can be neglected. How can the 
uncertainties be the same ? 

• The low energy scattering is expected according to the empirical fit of Ref. Q to lead to 0.862(12) fm, 
and one can obtain the high energy part from re-evaluation back from the average value 0.847(9) fm H]. 
It is close to 0.832(12) fm and that is it straight disagreement with the low energy result 1^. 

Concerning the second remark on the different uncertainties expected by us, we would like to say that the 
fitting has been done in Ref. for the isoscalar 



and isovector 



F[iq^)-Friq') 
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form factors, combined from the Dirac (i — 1) and Pauli (i = 2) form factors of the proton and neutron. 
We expect that the authors have obtained rather results for the fitting parameters but they are not ready to 
estimate properly the uncertainty for their specific function like the radii. 

Due to the different data from the low and high energy scattering we first discuss the different physics 
of the correction. 

• The two-loop QED correction are of different magnitude for the low and high momentum. 

o The low energy scattering is sensitive to them, 
o The high energy one is not. 

• The muonic and hadronic vacuum polarization are of the same relative 1% level for both of them. 

o The low energy data are sensitive. 

o The high energy ones are sensitive as well. 

• The analytic properties of the form factors can be of use. 

o For the Jow energy data evaluation they are not needed. It is a pure empirical fitting 
o The high energy data treatment in Ref. has been actually based on their using. 

• The normalization of the experimental data is the of one important points as we show below. 

o The low energy data are sensitive to that up to the level of 10~*. 

* In the Mainz-1980 work |^ it is included as a fixed parameter with the different values for 
different experiments. 

* In the Mainz-1996 paper ||] the normalization of the electric form factor is fixed absolutely to be 
equal to 1 at 5^ = . In this way one has to take care of the systematic errors due to correlation 
of the data from the same experiment. 

o In case of the high energy scattering, the data are sensitive to the normalization only on level up to 
10-2. 

• The magnetic form factor effects also have to be mentioned. 

o The low energy scattering is a mainly pure electric form factor scattering and magnetic effect are 
small. 

o The high energy data include such effects significantly. 

The short summary of the comparison 'low energy against the high one' can be described with following 
listing: 

• The low energy data are 

— sensitive to everything; 

+ needed no additional theory (e. g. like the dispersion approach); 
± only for Ge- 

• The high energy data are 

+ not so sensitive; 

— needed to be treated with using significantly a knowledge of the analytic properties; 
± both for the electric and magnetic form factors. 

The items labeled by '+' are advantages and the ones with ' — ' are disadvantages. 

Now after the comparison of the low energy and high energy data we would like to discuss some interpre- 
tation of the double Mainz difference. We have reported on some possible disagreement between the low 
energy and high energy data, but actually there are two options: 

• The disagreement can be actual an one and it has been discussed some possible corrections which are 
different for the high and low energy data. 
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• In the new Mainz work [|| they could obtain the different result from the one in Ref. for the low energy 
part of the data, because they could use another fitting. Some different details of their approaches have 
been also discussed above. 

Unfortunately in Ref. ||^ the authors have given no discussions which part of the data is mainly statistically 
responsible for the final result of the proton radius. If we suppose that the difference appeared from the low 
energy data, the main source of them is the experimental results of Ref. |^ . Maybe some clarifying remarks 
have been just done by Wong (ironically his work jssj was appeared in 1994 before the publication of the paper 
two years later). In the quoted work, only the Mainz data have been fitted according eq.(H), but in three 
different ways. We present all of them in Fig. 4. 

• The normalization constant is fixed to gq — 1.0014. The result is Rp — 0.863(9) fm close to Rp = 
0.862(12) fm from Mainz experimental paper ||]. The procedure is also close to the one of Ref. [^. The 
difference may come from including the data of the two other less precise experiments or from a different 
treatment of the systematic errors. 

• The normalization constant is fixed to uq = 1. The result is Rp = 0.849(9) fm close to Rp = 0.847(9) 
fm from the dispersion fitting work The same normalization was used there. 

• The normalization parameter qq has been treated as a free parameter to be found from the fit. The results 
are oq = 1.0028(22) and Rp = 0.877(24) fm. 

We think that only the last fitting can be correct, because the normalization constant ao cannot be measured 
with sufficient accuracy and has to be determined by fitting. It has to be mentioned here that in Ref. |^ a 
special run of measurements for the absolute gauging has been performed. The run included 5 points with 
the systematic errors equal to either 0.44% or 0.46% and the statistical errors in a range between 0.42% and 
1.28%. That means that the last fitting is actually in agreement with the experimental gauging within its 
0.7%-uncertainty and there is no reason to fix the normalization as it was done in the first two fits (and in Refs. 

The final uncertainty of the last result is significantly larger than of the two first. Such large changes of the 
uncertainty are due to two simple reasons. The fit is almost a straight line. To determine it one has to know 
a few points from the both edges of the data. For the straight line the value of the reference interval is quite 
important for the uncertainty When uq is fixed exactly one extra point at = is actually included. The 
reference interval is larger, because for the free normalization one has to use actual results of the cross-section 
measurement and in this case q±,„ = 0.14 fm"^. So we can see that for the free normalization the reference 
interval is now smaller and the number of the fit parameters is larger. 

5 Conclusions and discussions 

We can now summarize some results on the determination of the proton radius. To the values presented in 
Fig. 1 we have added in Fig. 5 five other values. First we have included without any discussion an old result 
with a more reasonable uncertainty (that is fitting in Ref. |^ of the Saskatoon data). Next we present in 
Fig. 5 the Wong fitting result with the free normalization approach. We also give three values based one the 
hydrogen Lamb shift study. One appears from using of the grand average value, and the others come from the 
absolute measurements. Due to the discrepancy mentioned above between the two Paris measurements we give 
the results for both 2s — 8s /d transition frequency values. 

5.1 What has been actually done 

• The Lamb shift value for the proton charge radius significantly depends on the data used. It is necessary 
to adjust all data. The uncertainty includes 0.010 fm (theory of the Is Lamb shift), 0.002 fm (theory of 
A(2)), 0.010 fm (expected experimental uncertainty). The final uncertainty is going to be 0.014 fm. 

• The lattice calculations cannot give any good result at least because only a few of the points for the form 
factor are found at moderate momenta. It is not enough for a successful fitting. We also have some doubts 
due to the uncertainty given in Refs. ||, |^. 

• The scattering data include an underestimation of the uncertainties. The best ones are from Mainz, and 
due to the Wong fitting with the free normalization it seems that the reasonable uncertainty has to be 
close to 0.024 fm (i. e. twice larger than in the original work |^). 
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• It is not clear what is the actual uncertainty of the radius from the Mainz dispersion fitting and which 
part of data have actually contributed to the result. 

5.2 What can be done 

• One of the most important values for atomic applications is — Maybe it is less dependent on 
the proton model. That problem has to be investigated. The value is to be helpful for extrapolating the 
low energy data, for extracting Rm from the Lamb shift measurement and for calculating of the Zemach 
correction (see Appendix for details). 

• Adjustment for the overall hydrogen data has to include also results for the proton radius. But the problem 
is: all results are correlated because obtained by fitting of the experimental data. The newer compilations 
of them are wider, but include also the older results. 

• Maybe a measurement of the muonic hydrogen Lamb shift can give a more safe value, but that is an 
open question. In case of the hydrogen-like ion of the helium-4, a high accuracy has been obtained in an 
experiment ||68|| , but the result has not been able to be reproduced by an other independent team [|69| . So 
the safety can be just a problem. 

• It is not actually clear what the true uncertainties of the Zemach correction with the Mainz dispersion 
fitting are. 
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A Hydrogen hyperfine structure and the proton structure 

The hyperfine splitting is much more sensitive to the proton structure than the Lamb shift. The main nuclear 
structure dependent contribution (so-called 'Zemach correction') is of the form 



Aiy{Ze'niach) — up 



2Zamp 



r.4 



1 + K 



1 



The Zemach correction includes integration of the form factors, but the most important part comes from the 
low momenta. That is because the high momenta asymptotic behavior of the integral is determined by '-1' and 
it is form factor independent. The low momenta asymptotics are defined by eqs. (|^) and (^ and the result for 
the Zemach integral is proportional to R^ + R^j. So, the result depends straightly on the both proton radius. 
The low momentum part of the integral comes from the momentum below 0.35 Gev ~ 3 fm~^ and one can see 
from Figs. 3 and 4, that in this region the approximation of eq.(^) can give an appropriate expression for at 
least the electric form factor. 

Some higher order structure dependent corrections were calculated by us in Ref. pst , and all old results 
were also reviewed there. That have allowed to give some prediction. According to a classical review by G. T. 
Bodwin and D. R. Yennie, 1988 [ |66| , the comparison of theory and experiment leads for the hydrogen hyperfine 
splitting to 

Vhfsiaxp) - Vhfs{theo) 



Vhfs{exp) 



0.48 ± 0.56 ppm 



Proton polarizability is not included in Vh fs (theo) and the difference above has to be interpreted as its contri- 
bution. The theoretical limitation for the proton polarizability contribution is pq 



Vhfsiexp) 



< A ppm. 
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That is obtained from a treatment of the inelastic scattering data. Our result for the comparison was 



Due to the status of the magnetic radius which is not good determined we need to say here that the uncertainty 
can be larger by factor like 3 from the pure extrapolation data in the experimental work Q . Using the dispersion 
fitting data Q the result has to be shifted slightly, but the uncertainty has to be approximately the same (0.9 
ppm), because Re = 0.847(9) fm and Rm = 0.853(9) fm and so the radii are equal one to the other within 
their uncertainty. According the Wong fit the result for the polarizability is rather close to 4(2) ppm within 
assumption, that the magnetic radius is equal to the electric one. 
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